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Abstract - We study the effects of time-varying environmental noise on nonequilibrium phase 
transitions in spreading and growth processes. Using the examples of the logistic evolution equa¬ 
tion as well as the contact process, we show that such temporal disorder gives rise to a distinct 
type of critical points at which the effective noise amplitude diverges on long time scales. This 
leads to enormous density fluctuations characterized by an infinitely broad probability distribu¬ 
tion at criticality. We develop a real-time renormalization-group theory that provides a general 
framework for the effects of temporal disorder on nonequilibrium processes. We also discuss how 
general this exotic critical behavior is, we illustrate the results by computer simulations, and we 
touch upon experimental applications of our theory. 


Introduction. — Systems far from thermodynamic 
equilibrium can undergo abrupt transitions between differ¬ 
ent nonequilibrium steady states. These nonequilibrium 
phase transitions feature collective behavior over long dis¬ 
tances and times just as thermodynamic equilibrium tran¬ 
sitions. Examples are found in the extinction dynamics of 
epidemics and bio-populations as well as in surface growth, 
turbulent flow, and catalytic reactions m- 

Nonequilibrium processes often occur in spatially inho¬ 
mogeneous systems and time-varying environments, i.e., 
in the presence of spatial and temporal disorder. Spatial 
disorder can have dramatic effects on nonequilibrium tran¬ 
sitions. For example, it destroys the ubiquitous directed 
percolation (DP) [5] universality class and produces an 
exotic infinite-randomness critical point mm. The asso¬ 
ciated Griffiths phases are dominated by rare fluctuations 
and feature anomalous power-law relaxation | Kill 171 . Sim¬ 
ilar behavior is found in percolating jlSlUfl] and quasiperi- 
odic [23] systems (for a review, see Ref. pH]). 

The effects of environmental noise, i.e., temporal disor¬ 
der, have attracted less attention. Kinzel [22] showed that 
temporal disorder destabilizes the DP transition because 
its correlation time exponent z/| = z/j \_z violates the gen¬ 
eralized Harris criterion zz|| > 2. Jensen [231121] employed 
series expansions and Monte-Carlo simulations to deter¬ 
mine the fate of the DP transition with temporal disorder 


and reported nonuniversal critical exponents. In an in¬ 
triguing paper, Vazquez et al. [25] demonstrated that rare 
noise fluctuations can lead to a temporal analog of the 
Griffiths phase, featuring an unusual power-law relation 
between lifetime and system size. Environmental noise 
has also been studied within space-independent (single¬ 
variable) models of biological population dynamics (see, 
e.g., Refs. pmiE])- Despite these efforts, a general frame¬ 
work for understanding nonequilibrium phase transitions 
in the presence of external noise is still lacking. 

In this Letter, we therefore study the effects of envi¬ 
ronmental noise on two prototypical models of spreading 
and growth processes: the space-independent (mean-field) 
logistic evolution equation and the spatially extended con¬ 
tact process. In both cases, we find highly unusual behav¬ 
ior close to the extinction transition. It is characterized by 
a diverging effective noise amplitude and enormous den¬ 
sity fluctuations on long time scales, motivating the name 
“infinite-noise critical point”. Infinite-noise critical behav¬ 
ior can be seen as counterpart of infinite-randomness crit¬ 
ical behavior in spatially disordered systems, but with ex¬ 
changed roles of space and time. 

In the remainder of this Letter, we first consider the 
logistic evolution equation where a clear picture of the 
infinite-noise physics emerges from an asymptotically ex¬ 
act random walk approach. To establish a framework for 
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the analysis of temporal disorder, we develop a real-time 
renormalization group (RG). We then apply the RG to our 
second model, the contact process, and find that the DP 
critical behavior of the pure problem gets replaced by an 
infinite-noise critical point as well. We confirm key find¬ 
ings by computer simulations, and we discuss applications. 

Logistic equation. — The logistic evolution equation 


describes a variety of growth processes in nature, with 
applications as diverse as biological population dynamics 
[29] , catalytic chemical reactions [30] , and even linguistics 
[31] . In the context of epidemic spreading, p is the density 
of infected individuals, p is the rate at which sick indi¬ 
viduals heal while A corresponds to the rate at which a 
healthy individual is infected by sick ones. If A and p are 
time-independent, the behavior is easily understood; the 
epidemic dies out exponentially for A < p and survives for 
infinite time for A > p. At the critical point, A = p, the 
density decays to zero, but only as a power law, p ~ 1/t. 
Environmental noise, i.e., temporal disorder, can be in¬ 
troduced by making the healing and infection rates time 
dependent. For definiteness, we consider the rates to be 
piecewise constant, p(t) = p n and \{t) = X n , over time 
intervals A t n . The p n and \ n are independently drawn 
from probability distributions W fl (p) and WA(A). 

If the healing and infection rates are time-independent, 
eq. (H]) can be solved in closed form. Using this solution 
within each time interval of constant rates yields a linear 
recurrence for the inverse density, 

Pn+1 ~ a nPn + c n • (2) 

p n is the density at the start of time interval n. The 
multipliers a n = exp[(/x„ — X n )At n ] reflect the exponential 
growth or decay due to the linear term in eq. (jX)). The 
constants c n = {a n —l)X n /(p n — X n ) are only important for 
large densities; they limit the growth and prevent p n > 1. 
The time evolution is thus a random sequence of decay 
and spreading segments during which the density p either 
decreases or increases, depending on the balance between 
p{t) and A (t). This is illustrated in Fig. [T] Neglecting the 
c n in the recurrence (O for the moment, we see that the 
logarithm of the density, x n = — In p n , performs a simple 
random walk, x n+ i = x n + In a n . The effect of the c n can 
be approximated by a reflecting boundary for this random 
walk at x = 0. It limits the x n to positive values. 

The reflected random-walk theory allows us to find the 
time evolution of the full probability distribution of the 
density rather than just its average. If the evolution starts 
from a fully infected system, po = 1, the probability dis¬ 
tribution P n (x) after a large number n of time intervals is 
given by (for x > 0) 


P n (x) = 
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Fig. 1: (color online) Density p{t ) of individual noise realiza¬ 
tions, plotted as x = — In p vs. t. The logistic evolution results 
are obtained from eq. (1) with A = 1, At = 1 and binary disor¬ 
der W p (p) = pS(p-ph) + {l— p)S(p-pi) with pn = 2 ,pi = 0.5 
and p — p c = 1/3. The 2d contact process data stem from a 
Monte-Carlo simulation of 1000 2 sites using p = 1, At = 2 and 
Wx(A) = pS(X — Ah) + (1— p)S(X—Xi) with Xh = 3.65, A; = 0.365 
and p = 0.8. Inset: Schematic of the RG. In each step, the seg¬ 
ment with the smallest change of x = — In p is eliminated. 


Here, $ is the cumulative normal distribution, v = (In a) 
yields the drift of P„(x), and a 2 = (In 2 a) — (In a) 2 gives its 
width. [(...) denotes the average over W p (p) and Wa(A).] 
The solution © can be verified by inserting P n (x) into 
the drift-diffusion equation d n P = (a 2 /2)d 2 P—vd x P with 
flux-free boundary condition vP—(a 2 /2)d x P = 0 at x = 0. 

By inspecting the distribution P n {x ), we can identify 
three different regimes: v < 0 corresponds to the active 
phase in which the epidemic survives for infinite time, v > 
0 is the inactive phase in which the epidemic dies out, and 
v = 0 is the critical point separating the two. v thus 
serves as a measure for the distance from criticality. We 
now discuss the regimes in more detail. 

In the active phase, v < 0, P n (x) approaches the sta¬ 
tionary distribution P(x) = {2\v\/a 2 ) exp[— 2\v\x/<j 2 ] in 
the long-time limit. The distribution P p (p) = p~ 1+1 / K /n 
of the density itself is highly singular and characterized 
by a non-universal Griffiths exponent k = cr 2 /(2\v\) that 
diverges at criticality. The average stationary density is 
given by (p) = (exp(— x)) = 1/(1 + re) while the typical 
density reads p typ = exp(—{x)) = exp (—re). Close to crit¬ 
icality, the average density (which is dominated by rare 
events) is much larger than the typical one. 

At the critical point, v = 0, the distribu¬ 

tion P n {x ) simplifies to a half-Gaussian, P n (x) = 
2(27r<T 2 n) -1 / 2 exp[— x 2 / {2a 2 n)\, that broadens without 
limit with increasing time. This clearly illustrates the 
notion of infinite-noise criticality. For long times, the 
average density thus decays as (p n ) = (exp(— x n )) = 
2(27r<r 2 n) -1 / 2 ~ t -1 / 2 . In contrast, the typical 
density decays much faster, p typ = exp(— (x n )) = 

exp[— (2<7 2 n/7r) 1/,2 ] ; implying In p typ ~ — t 1 / 2 . 
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The inactive phase, v > 0, is more conventional as the 
entire distribution P n (x) moves to larger x with increas¬ 
ing time. Thus, the density of almost all noise realiza¬ 
tions rapidly vanishes such that both average and typ¬ 
ical densities decay exponentially with time. Specifically 
p typ ~ exp(— vn) while (p n ) ~ exp[— nv 2 /{2a 2 )] for v < a 2 
and (p n ) ~ exp[— n{v — a 2 / 2)] for v > a 2 . The correlation 
time is given by the time when the off-critical P n {x) 
starts deviating significantly from its critical counterpart. 
This yields = {a 2 /v 2 )At. 

We can cast these results in the language of critical phe¬ 
nomena by comparing our results with the definitions of 
the critical exponents, (p) ~ t~ s at criticality, (p) ~ \v\& 
in the active phase, and £t ~ M _!/|1 • The exponent values 

8 = 1/2, 0 = 1, u ||=2. (4) 

differ from the “clean” ones (6 = 1, 0 = 1, v\\ =1 [2]) 
but fulfill the scaling relation 5 = 0/v\\, and vu saturates 
Kinzel’s bound z^j > 2 [22] , What about the correlation 
length? The logistic equation does not contain any notion 
of space. However, if it describes a d-dimensional system 
of N individuals, we can introduce a length scale L ~ 
N 1 / d . It is clear that the critical dynamics changes when 
the typical density p typ = exp[—(2cr 2 n/7r) 1 / 2 ] becomes of 
order of 1/N. This suggests an exponential dependence 
between correlation length £ and correlation time £ 4 , 

ln£ ~ £" with oj = 1/2. (5) 

The dynamical exponent z is thus formally zero, and the 
correlation length exponent = v\\/z is infinite. This 
highly unusual dynamical scaling is analogous to the acti¬ 
vated scaling at infinite-randomness critical points in spa¬ 
tially disordered systems mm , but with the roles played 
by space and time exchanged. 

To illustrate and verify the reflected-random-walk the¬ 
ory, we have solved the logistic equation (ID numerically. 
Figure [2] shows the resulting P n (x) and ( p(t )). All data 
confirm our predictions. We have also verified the station¬ 
ary distribution P(x) = (2|u|/ct 2 ) exp[— 2|u|:r/cr 2 ] in the 
active phase close to criticality. 

The knowledge of the full density distribution allows us 
to obtain many additional results. The typical lifetime tn 
of a finite system of N sites in the active phase can be 
estimated from the probability that the density is below 
1/N. This yields r^ 1 ~ f^ N dpP p {p) ~ -ZV _1//re . At crit¬ 
icality, the condition p typ (rAr) = 1/N gives tn ~ ln 2 (iV), 
and in the inactive phase we find tn ~ ln(7V). All these 
lifetimes agree with the Langevin results of Ref. EH] . 

Real-time renormalization group. — To establish 
a general framework for treating the temporal disorder, 
we now develop a real-time RG, in analogy to the strong- 
disorder RG [34] for spatially disordered systems. We first 
combine consecutive time intervals with p > A into a single 
interval of length At up (and consecutive intervals with p < 
A into a single interval of length Af dn ). This leaves us 



Fig. 2: (color online) Numerical solution of the logistic equation 
© for A = 1, At = 1 using 10 5 disorder realizations drawn from 
W(p) = pS(p — ph) + {l—p)S(p — pi) with pn = 2 ,pi = 0.5. The 
transition is tuned by p, the critical point is at p c = 1/3. Main 
panel: Probability distribution P n (x) at criticality for several 
times t. The solid lines are the predictions of the random-walk 
theory, without adjustable parameters. Inset: Average density 
(p) vs. time t for several p. The dashed line is a fit to the 
prediction (p) ~ t -1 ? 2 . 

with a zig-zag curve of alternating upward and downward 
segments for x = — In p, as illustrated in the numerical 
results of Fig. [T] In each segment, p -1 follows the linear 
recurrence p~ 1 (t + At) = ap _1 (t) + c with the multipliers 
of the upward segments fulfilling a up > 1 and those of the 
downward segments a dn < 1. 

The idea of the RG consists in iteratively eliminat¬ 
ing the smallest upward and downward segments, coarse- 
graining time. The distributions of a up and a dn change 
during this process, and their asymptotic behavior gov¬ 
erns the long-time physics. In each RG step, we find the 
multiplier a closest to unity. This defines the RG scale 
Q = min(a“ p , l/a dn ). We then replace this segment and 
its two neighbors by a single segment, as shown in the 
inset of Fig. |T] The time evolution of the density in this 
renormalized segment follows from combining the three 
original density recurrences. It takes the same linear form 
p _1 (t + At) = ap~ l (t) + c but with renormalized coeffi¬ 
cients a and c. The renormalized multipliers read 


d up = (6) 

l/a dn = (l/a dn )(l/a d " 1 )/H, (7) 

for the decimation of downward segment a dn and upward 
segment a“ p , respectively. The time intervals combine as 

At up = Af“ p + Af dn + At“ p x , (8) 

At dn = Atf^ + Af“ p + At dn . (9) 


The expressions from |[6|) to |9|) are exact. Moreover, they 
are equivalent to the RG for the spatially disordered quan¬ 
tum Ising chain mm and can be solved in the same way. 

Following Fisher mm , we define logarithmic variables 
r = In fl, 0 = In a up — F, and £ = — In a dn — F. Upon iter- 
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ating the RG steps, T increases, and the probability dis¬ 
tributions of f3 and 75(/3;T) and 7FC>r) change. They 
fulfill the integro-differential (flow) equations 

If = ^ + (Fo-'Po)K + v 0 (n&nj . (io) 

dV dV ( C \ 

OF = -g^ + ^o-^V + TZoiv^V) . ( 11 ) 

where TZo = 75.(0;T) and Vo = V(0 ;T). The symbol V ® 
V = fg V((')V(( — C)d(' denotes the convolution. 

The complete solution of the flow equations for 75(/3; T) 
and V((’,T) is quite elaborate [32|j33j|35]. However, the 
physically relevant solutions can be obtained from the 
ansatz 

75(/3; T) = 7Z 0 e~ n °P , V{(; T) = Voe~ v ^ . (12) 

Inserting this ansatz into the flow equations m and m 
gives the ordinary differential equations 

cLIZo/dr = — IZqVq , dVo/dV = — TZqVq ( 13 ) 

which can be solved by elementary means. 

This leads to three types of solutions for T —oo. At 
criticality, (£) = (/3) (equivalent to v = (In a) = 0), 
the distributions of the upward and downward multipli¬ 
ers are identical, T^CiT) = exp(—£/r)/T and 75(/3;T) = 
exp(—/3/T)/T. They broaden without limit with the RG 
scale T, i.e., the effective temporal disorder diverges for 
long times, in agreement with the notion of “infinite- 
noise” criticality. The typical length of a renormalized 
time interval At behaves as T 2 . In the inactive phase, 
v = (/?) — (C) > 0, the distribution of the upward 
steps becomes very broad, 75.(/3; T) = IZo exp(— 1Zq(3) with 
TZo = exp(— T/k)/k while the distribution of the down¬ 
ward steps does not broaden, V(Ci T) = exp(— (/k)/k with 
a constant k ~ 1/v. In the active phase, v < 0, the ex¬ 
pressions for ^(Cjr) and 75(/3;T) are exchanged. In both 
phases, At scales as exp(r/«). 

Many important results follow from the RG. To estimate 
the lifetime of a finite-size system of N sites, we run the 
RG until the typical upward multiplier a up = exp(T + ft) 
reaches N. The lifetime tjv is given by the typical time 
interval at that RG scale. Using the above solutions, we 
obtain tn ~ N 1 /' t in the active phase, tn ~ (In A) 2 at 
criticality, and tn ~ In A in the inactive phase, as in the 
random walk approach. Analogously, the width of the 
density distribution P{x) at criticality is governed by the 
upward multipliers a up . As the typical value of lna up = 
T + /3 is 2T, we obtain a width of Ax ~ T ~ t 1 / 2 in 
agreement with Q 0. 

Contact process. — So far we have analyzed the lo¬ 
gistic equation (JX|) which contains density fluctuations in 

1 The full P(x) can be found from the joined distribution of a and 

c of the last segment at the end of the decimations. 


time but not in space. What changes for finite-dimensional 
systems that fluctuate in space and time? To answer this 
question, we turn to the contact process [36] . It is defined 
on a d-dimensional lattice. Each site can either be infected 
or healthy. During the time evolution, each infected site 
heals at rate p while healthy sites becomes infected at 
rate Xn/(2d) where n is the number of infected nearest 
neighbor sites. (Note that the logistic equation dT]) can be 
understood as the mean-field limit of the contact process.) 
If /i and A are uniform in space and time-independent, the 
nonequilibrium transition between the active and inactive 
phases is in the ubiquitous DP universality class 0 . 

Temporal disorder is again introduced by making the 
healing and infection rates random functions of time; we 
assume piecewise constant rates, n(t) = p n and A(f) = A„, 
over time intervals A t n . According to Kinzel’s [22] gen¬ 
eralization of the Harris criterion, such temporal disorder 
destabilizes the DP critical behavior in all dimensions be¬ 
cause the clean DP correlation time exponent z/|| = v±z 
violates the inequality i/|| > 2. Temporal disorder is thus 
a relevant perturbation (having positive scale dimension) 
at the clean DP critical point. 

To resolve the fate of the phase transition in the pres¬ 
ence of temporal disorder, we generalize our real-time RG 
method. As in the logistic equation, the time evolution of 
the contact process is a sequence of spreading and decay 
segments, depending on the values of n n and A„. If the 
temporal disorder is strong, the system is far away from 
criticality in each individual segment. This allows us to 
neglect spatial fluctuations and to formulate the theory 
in terms of the density p(t) only. Wc will discuss the va¬ 
lidity of this approximation in the conclusions. During 
decay segments, the density decreases exponentially with 
time, as in the (mean-field) logistic equation, because each 
site can heal independently. In contrast, the character of 
the spreading segments changes because in a finite-d sys¬ 
tem with short-range couplings, the infection can at best 
spread ballistically, p(t ) ps po(l + bt) d , rather than ex¬ 
ponentially as in the logistic equation. This difference is 
clearly visible in Fig. [T] (This also implies that x = — In p 
does not undergo a simple random walk, consecutive steps 
are rather correlated in a nontrivial manner.) 

The RG must thus be modified. If a downward segment 
(density increase) is decimated, the multipliers a up still 
renormalize multiplicatively according to |6]). If an up¬ 
ward segment is decimated, we need to combine the two 
neighboring downward (ballistic spreading) segments dur¬ 
ing which the radii of active clusters grow linearly with 
time. For strong disorder, l/af n , 1/afi^ H, we can 
therefore estimate the renormalized multiplier as 1/a dn = 
(1 + bAt dn ) d w (bAt dn ) d = (biAt dn + 6 i _iA^“ 1 ) d ps 
[{l/a dn ) 1 / d + (1 /a d * 1 ) 1 / d ] d . We thus arrive at an addi¬ 
tive rather than a multiplicative renormalization 

(l/d dn ) 1/d = (l/a dn ) 1/d + (l/a d ^ 1 ) 1/d -H 1 / d (14) 

Here, the last term contains the subleading contribution 
of the decimated upward segment which we have mainly 
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added to ensure the correct behavior in the atypical case 
l/a dn = = Q, The RG defined by recursions © 

and (fTdl) is equivalent to that of spatially disordered quan¬ 
tum systems with super-Ohmic dissipation [571 or with 
long-range interactions [38] . @ 

To analyze the RG, we define reduced variables 0 = 
lna up — T and £ = d[(fla dn ) -1 / d — 1], In terms of these 
variables, the flow equation for the distribution TZ(j3\ T) is 
still given by (flTil) while V(C;T) fulfills 


dV 

dT 


= 1 + 


&P 

d( 


+ ('Po — n 0 h —-) v +izq (v ® v 


(15) 

In the mean-field limit d —> oo, we recover the flow equa¬ 
tion (ED. as expected. Inserting the exponential ansatz 
tm> into the flow equations m and m yields 


dllo/dT = -K 0 Vo , dV 0 /dT = (1/d - TZ 0 )P 0 . (16) 


These equations take the famous Kosterlitz-Thouless form 
[39] for all finite d. Details of their solution will be pub¬ 
lished elsewhere. Here, we just summarize the results. The 
critical fixed point is located at Vq = 0,1Zq = 1/d. This 
implies that the distribution of the downward multipli¬ 
ers l/a dn becomes infinitely broad while the distribution 
upward multipliers a up retains a finite width for T —> oo. 
The typical length of a renormalized time interval behaves 
as At ~ = Vl l / d . 

Many important results follow. The width of the (loga¬ 
rithmic) density distribution P(x) at criticality is governed 
by the typical value of lna up . This gives Ax ~ lna up « 
T = lnf2 ~ lnt. Thus, P(x) still broadens without limit 
with increasing time, but its width increases only loga¬ 
rithmically, in contrast to the power-law increase in the 
case of the logistic evolution equation. The average den¬ 
sity thus decays as (p(t)) ~ (In t)~ s with <5 = 1 (while the 
usual decay exponent <5 formally vanishes). The correla¬ 
tion time can be obtained from analyzing the off-critical 
solutions of (fTBl) . As usual for Kosterlitz-Thouless flows, 
it depends exponentially on the distance from criticality 
r via ln£t ~ with P|| =1/2 (the usual correla¬ 

tion time exponent v\\ is infinite, fulfilling Kinzel’s bound 
v\\ > 2). Correlation length and time are proportional to 
each other, implying a dynamical exponent z = 1. Fi¬ 
nally, the stationary density in the active phase varies as 
\r\P with /3 = 1/2. The critical exponents 

S = l, 0=1/2, P|,=l/2 (17) 

fulfill the scaling relation 6 = fd/v\\. In the active phase, 
the lifetime r jv of a finite-size sample shows a Griffiths 
singularity, i.e., it increases as a nonuniversal power law 
Tat ~ IV 1 / K . However, in contrast to the logistic evolu¬ 
tion equation, the Griffiths exponent k does not diverge 
at criticality but saturates at k c = d. 

2 If the density increase in the spreading segments were slower 
than ballistic but still followed a power in t, the recursion for l/a dn 
would take the additive form JED with 1/d replaced by a different 
exponent. The resulting critical behavior would not change E2- 



Fig. 3: (color online) Monte-Carlo simulations of the 2d contact 
process, starting from p = 1 (up to 3200 2 sites, p = l,At = 
2, W(A) = 0.8 5(\ - A h ) + 0.2 J(A - Ah/10), 20000 disorder 
realizations). Main: P(x) at criticality for different times t, 
scaled such that the curves coincide. Inset (a): Scale factor 
ft vs. In L. confirming the logarithmic broadening of P(x) with 
time. Inset (b): (p(t)) for varying infection rate A. The critical 
curve, Xh = 3.65, follows the predicted (p) ~ (lnt)^ 1 over 
almost three orders of magnitude in t. 


To test these predictions, we have performed Monte- 
Carlo simulations of the 2d contact process with temporal 
disorder by methods analogous to Ref. m- The time 
evolutions of P{x) and (p(t)) are shown in Fig. [3] All 
data are in excellent agreement with the RG predictions. 
We have also verified lyv ~ iV 1 / 2 at criticality. 

Conclusions. — In conclusion, random environmen¬ 
tal noise (i.e., temporal disorder) leads to exotic “infinite- 
noise” critical points at which the density distributions be¬ 
come infinitely broad (even on a logarithmic scale) both 
for the (mean-field) logistic evolution model and for the 
finite-d contact process. The infinite-noise critical behav¬ 
ior of the logistic equation is asymptotically exact because, 
at criticality, the typical density vanishes with increas¬ 
ing time, justifying the reflected-random-walk approach. 
Moreover, the real-time RG suggests that this critical be¬ 
havior holds for an entire class of single-variable growth 
models with temporal disorder in the linear growth rate 
and some nonlinearity preventing population explosion. 
Such models occur in population dynamics, chemical ki¬ 
netics, economics, and other fields. 

The generalization of the RG to finite dimensions in¬ 
cludes spatial fluctuations only approximately (via the 
modification of the density increase in spreading seg¬ 
ments). The approximation is expected to be good if the 
temporal disorder is so strong that individual spreading 
and decay segments are far from criticality. As the density 
distribution broadens without limit at criticality, this con¬ 
dition appears to be fulfilled self-consistently. Moreover, 
according to Kinzel’s [22] generalization of the Harris cri¬ 
terion, even weakly disordered systems should flow to this 
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regime. These heuristic arguments suggest that our the¬ 
ory is asymptotically stable for all bare disorder strengths. 
However, two other scenarios cannot be excluded: (i) our 
theory might be asymptotically stable only for sufficiently 
strong bare disorder, or (ii) our theory might hold in a 
transient time-regime (whose width diverges with disorder 
strength) before spatial fluctuations become important. 
Discriminating between these scenarios requires more so¬ 
phisticated methods. 

Our theory for the contact process predicts universal 
critical behavior while Jensen (23)124) reported nonuni- 
versal exponents in Id. (Jensen’s values do not always 
respect the bound U|| > 2, though.) This could mean 
that our theory does not apply to Id, or Jensen’s data 
could be preasymptotic, caused by long transients before 
the infinite-noise behavior is reached. We have performed 
Monte-Carlo simulations of the Id contact process with 
strong temporal disorder; preliminary results do show sig¬ 
natures of the strong-noise physics reported here. 

In recent years, absorbing state phase transitions have 
been observed in turbulent liquid crystals [40] , driven sus¬ 
pensions [4UI42] . superconducting vortex dynamics [43], 
as well as in bacteria colony biofilms mm- Studying ex¬ 
ternal noise at these transitions can provide experimental 
tests of our theory. Moreover, the effects of noise on the 
extinction of a biological population or an entire species 
due to environmental changes are attracting considerable 
attention (see, e.g., Ref. [21]). Experimentally, these ques¬ 
tions could be studied, e.g., by growing bacteria or yeast 
populations in fluctuating external conditions. 
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